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ene DUC TION 


foe PURPOSE AND HYPOTHESIS 

The purpose of this study is to test the applicability 
of the Garwood (1977) one-dimensional bulk mixed layer 
model for hindcasting the cold ocean thermal anomalies 
found in the Anomaly Dynamics Study (ADS) domain. Elsberry 
(1983) describes numerous cold anomalies and selects 
several anomalies that are particularly suitable for 
testing ocean prediction models. Two of these cold 
anomalies were chosen for this study. The ability to pre- 
dict departures from climatology provides a useful test of 
a model's capabilities relative to either a persistence or 
a climatological (zero anomaly) forecast. The study of 
anomalies also lends insight into the large-scale variability 
of the ocean thermal structure. A second purpose of this 
study is to demonstrate the usefulness of the corrections 
derived by Elsberry et al. (1982) to Fleet Numerical 
Oceanography Center (FNOC) surface heat flux estimates in 
North Pacific Ocean predictions. ThesSe corrections are 
necessary to offset a systematic bias in the FNOC heat 
fluxes which was discovered in the prediction experiments 
Memes berry et al. (1979), Budd (1980) and Steiner (1981). 

The seasonal variation in heat content of the upper 


ocean, away from the major current systems, is primarily 
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determined by the net heat flux across the ocean-atmosphere 
interface. This heat is distributed in the vertical almost 
exclusively by turbulent mixing. During and after the 
formation of the seasonal thermocline, the gradual increase 
in the solar flux tends to be offset by upward surface heat 
fluxes and increases in entrainment mixing associated with 
strong atmospheric storms. That is, a decrease in sea- 
surface temperature is found during periods of higher wind 
speeds (Elsberry and Garwood, 1978). These higher wind 
events occur less frequently during the spring and summer. 
The balance between periods of net warming and net cooling 
is in favor of increasing sea-surface temperature from the 
spring transition until early autumn, After this time, 
the net surface cooling and entrainment mixing associated 
with strong storms reverses the balance, and the sea-surface 
temperature decreases. Superposed on these seasonal trends 
are periods with above or below normal temperatures. 
Anomalous sea-surface temperatures can be caused by 
anomalous solar radiation and surface heat fluxes or by 
anomalous entrainment heat flux at the mixed layer base 
generated by wind stirring and convective overturning 
(Elsberry and Garwood, 1978) and by horizontal advection, 
Ekman pumping and upwelling. The basic hypothesis of this 
study is that near-surface cold temperature anomalies in 
the North Pacific Ocean during 1976-1978 were primarily 


generated by local vertical mixing processes. 
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PeeeeolUDY DESCRIPTION 

An objective of the North Pacific Ocean Experiment 
(NORPAX) was the study of the interaction between large- 
scale ocean temperature anomalies and weather. One observa- 
tional component of NORPAX was a ship-of-opportunity 
expendable bathythermograph (XBT) program (TRANSPAC) 
designed to observe the ocean thermal structure on space 
scales of thousands of kilometers (White and Bernstein, 
1979). The TRANSPAC analyses provide the initial and 
verifying temperature profiles that are required to validate 
ocean prediction models (Elsberry and Garwood, 1980). 
Because all of the TRANSPAC observations made within a 
particular month are used in the objective analysis, this 
defines the time and space scales for the initialization and 
verification of the model. 

The NORPAX ADS area (Fig. 1) is the oceanic region 
examined in this study. The ADS region is bounded by 30°N- 
50°N and 130°W-160°E and encompasses midlatitude regions 
that contain large-scale thermal anomalies. This region 
also has strong atmospheric variability. In this study, two 
cold anomalies are examined. The first is a large-scale, 
long-duration event that began in the fall of 1977, and 
extended into the spring of 1978 in the vicinity of 36°N, 
150°W. The second anomaly is a shallow, large-scale, rapid 
transition event occurring in the summer of 1976 in the 


vicinity of 40°N, 165°E. These two anomalies provide an 
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opportunity to test the Garwood model and the vertical 


mixing hypothesis. 


feeb’ INITIONS 

Several terms that will be used throughout this thesis 
are defined and discussed in this section. 

iw opring, Transition 

The transition from a winter mixed layer regime to 
a summer regime occurs during the spring, when the occur- 
rence and magnitude of storms (and thus high wind speeds) 
are diminishing and the net daily insolation is increasing. 
The increase of absorption of solar radiation in the near- 
surface layers promotes stability. 

Tully and Giovando (1963) noted that the spring 
transition appeared to be rapid. Elsberry and Garwood (1978) 
combined an oceanic model simulation and observations to 
demonstrate that the spring transition can occur during a 
Single daytime heating cycle. Budd (1980) reported that 
the spring transition usually occurred within a 36-hour 
period; however, the time period for transition was longer 
for steady wind conditions. The key synoptic feature ini- 
tiating the transition is an extended period of light winds 
cOinciding with a period of net downward heat flux (Elsberry 
and Raney, 1978). A layer of warmer (less dense) water 
near the surface is established during the diurnal heating 
cycle. If the mechanical generation of turbulent kinetic 


energy is sufficiently small, the layer will remain intact 


18 





fjmeeugn the night. A repetition of this cycle through the 
next several days can lead to the formation of the seasonal 
thermocline. 

A set of criteria based on the mixed layer depth 
and temperature was used to determine when the ocean 
boundary layer changed from a winter to a summer regime 
(Similar to Budd, 1980). The transition period was defined 
as the first period of sustained shallow mixed layer depths 
(<40 m) that followed a period of greater than 80 m depths. 
Occasionally the predicted and/or analyzed depths may later 
exceed 80 m for a short period a week or more after the 
establishment of the stable layer. Consequently, the 
transition period was specified as that period which also 
coincided with a significant increase in mixed layer tempera- 
ture. As the TRANSPAC temperature analyses were available 
monthly, only an approximate time frame can be given for the 
spring transition. 

2. Mixed Layer Depth 

There 1S no consensus method for determining the 
mixed layer depth of a vertical temperature profile. Various 
definitions have been applied in different studies depending 
on the requirements (operational or research) and profile 
resolution. In this study, the mixed layer depth has been 
defined to be that depth at which the vertical temperature 


gradient exceeds 0.10°C/5m. For regions with a very slight 
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negative gradient above a marked thermocline and for which 
the criteria of 0.1°C decrease in 5 m was too stringent, 
the mixed layer depth was manually determined from the 
vertical temperature profile. This was done to reduce 

any bias in mixed layer depth due to definition as found 
by Steiner (1981). 

In the Garwood model, a minimum mixed layer depth 
of 5 mwas set. This iS consistent with the analyses, as 
the wake turbulence of a ship from which an XBT was launched 
would destroy an shallow (<5 m) surface layer(s), and thus 
these layers would not be present in the analyses. The 
maximum model mixed layer depth was set at 190 m to prevent 


the model vertical profiles from becoming isothermal. 
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II. DATA SOURCES AND PROCEDURES 


A. DATA SOURCES 

This study is based on data between January 1976 and 
December 1979 from two separate sources. The FNOC atmo- 
spheric forcing and TRANSPAC temperature analyses are over 
the ADS domain bounded by 30°N, 50°N, 130°W and 160°E. The 
Seta resolution is 2° latitude by 5° longitude. 

Atmospheric forcing data (wind speed and direction, 
solar heat flux and total surface heat Flux) were obtained 
from the archived FNOC northern hemispheric atmospheric 
analyses and short-term predictions. The wind components 
were available at six-hour intervals, and the solar and 
total surface heat (latent plus sensible plus back radia- 
tion minus Solar) flux values at 12-hour intervals. The 
Garwood model requires atmospheric forcing fields on hourly 
time scales to resolve the oceanic response to the diurnal 
heating cycle. Garwood (1977) found that this diurnal 
component can modulate the seasonal trend. A complete 
description of the procedures and programs for editing and 
interpolating the forcing fields to one-hour intervals is 
tmecallacher (1979). 

The TRANSPAC data sources and analysis procedures are 
described in White and Bernstein (1979). XBT observations 


were made by ship-of-opportunity personnel along ship 
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tracks between Japan and the west coast of the United 
States. Due to ships tending to avoid rough weather there 
is a seasonal and a "fair weather" bias in the data. During 
the summer, the ship tracks are farther north and there is 
good data coverage over the domain. However, the northern 
part of the domain is not well sampled during the winter. 
It was found that the northwest and the southwest corners 
of the grid were not well sampled (Elsberry, 1983). Conse- 
quently, these areas are excluded in this study and will 
appear as cross-hatched areas on all appropriate figures. 
The XBT profiles have been objectively analyzed at 
@emchs of 0, 30, 60, 90, 120, 150 and 200 m. These 
objectively analyzed fields produced by White and Bernstein 
(1979) are not verticaliy coupled. Since the analyzed 
temperature profiles are available monthly, they are taken 
to be representative of the "observed" temperature profile 
@mmene 15th of the month. To assure a continuous record, 
wherever minor data gaps existed along the edge or in the 
interior of the ADS grid, time and space interpolations 
were made at 0, 60, 120 and 200 m (Elsberry, 1983). Since 
data gaps were not removed in the analysis at the other 
analyses levels, only the analyses at 0, 60 and 120 m were 
used to describe the horizontal temperature anomaly 
Structure. The horizontally edited values at the above 
levels, plus objectively analyzed values (if present) at 


30, 90 and 150 m, were then vertically interpolated to a 
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5 m Spacing between 0 and 200 m. These gridded monthly 
oceanic temperature fields were used to initialize the model 
Simulations and to verify the model output. 

The following comments should be kept in mind during 
subsequent discussions of the analyses over the entire 
domain. Due to the sparse ocean data, differences between 
the Garwood model predictions and the TRANSPAC observations 
of less than 0.5°C are probably not significant (Budd, 1980). 
The baseline study conducted by White and Bernstein (1979) 
indicated that the region west of 175°W was dominated by 
strong mesoscale (309 km) baroclinic eddies or waves, whereas 
east of 175°W the larger scale variability dominated. The 
temporal and spatial data coverage and the 2° latitude by 
5° longitude grid interval are inadequate to resolve eddies 
mm the Ses Hot west of the dateline and in the California 
current. Furthermore, the analyses in regions north of 
45°N and south of 35°N are less reliable due to the seasonal 


variation in data coverage (Elsberry, 1983). 


B. MEAN AND ANOMALY FrIELDS 

The monthly mean temperature (climatological) fields 
used in this study are an average of the TRANSPAC monthly 
objectively analyzed fields over the four years, 1976-1979. 
These mean temperature fields are based only on the TRANSPAC 
vertical temperature profiles, and do not include any other 


source of temperature observations. The advantage of the 
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four-year mean is that the analysis method is consistent 

in both the mean and the anomaly fields (Elsberry, 1983). 
The disadvantage is that the anomaly fields will differ 
from those generated using the long-term NORPAX mean fields 
(Barnett and Ott, 1976). This is particularly true because 
two of the four winter periods (1976 and 1977) include 
extreme conditions in the Central Pacific (White and 
Bernstein, 1979; Elsberry et al., 1979; White et al., 
meu; Haney, 1980; Budd, 1980). Thus cold anomalies during 
these winters have a smaller magnitude in this study than 
in the NORPAX fields. Elsberry (1983) found that there is 
an overall trend toward a warm bias in the four-year surface 
TRANSPAC mean fields relative to the long-term Barnett and 
@ee (1976) fields. 

The anomaly fields were formed by subtracting the four- 
year mean (climatclogical) fields from the analyzed fields 
at 0, 60 and 120 m. The resulting horizontal anomaly fields 
are taken to be the anomalous conditions on the 15th of the 


month. 


fee eRe DICTION MODEL 

The Garwood (1977) ocean mixed layer model is a second 
order closure bulk model. Since the model is one-dimensional, 
no advection or Ekman dynamics are included. A more complete 
description of the Garwood model is contained in the 
Appendix. 
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The Garwood model is initialized with the initial 
(TRANSPAC objectively analyzed) temperature profile fields 
for the 15th of the month. The model is used to predict 
the evolution of the oceanic thermal structure over the 
ADS area due to local surface forcing. The model predicts 
vertical temperature profiles in 5 m increments to 200 m 
over the entire grid every hour. For this study, the 
model output consists of the vertical temperature profiles 
averaged over 5 days centered on the 15th of each month, 
and the mixed layer temperature at the time of maximum 
mixed layer depth for each day. 

The model anomaly fields were formed by subtracting 
the four-year mean (climatological) fields from the model 
predicted fields at 0, 60 and 120 m. The resulting hori- 
zontal model anomaly fields are taken to be representative 


of the conditions on the 15th of the month. 
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III. LARGE-SCALE LONG-DURATION CASE 


A. COLD ANOMALY 12 DESCRIPTION 

The cold anomaly (CA) described below was designated 
and described by Elsberry (1983) as CA 12. CA 12 is a 
large-scale, long-duration event which persisted for 
approximately eight months and encompassed a vast region 
of the ADS grid. 

The entire domain had above average temperatures during 
October 1977 (Fig. 2A). For orientation purposes, the 
ultimate location of maximum intensity or central point 
(36°N, 150°W) for CA 12 is marked as a dot on the appli- 
cable figures. Starting in November 1977 (Fig. 2B), CA 12 
began to develop at the surface in the northeast region. 
The temperature decreases in this region are noteworthy 
relative to the climatological cooling during this period. 
In the northeast region, which changed from a warm to cold 
anomaly between October and November 1977, the surface 
temperature decrease was greater than 1.5 times the 
climatological decrease (Fig. 3). The cold anomaly con- 
tinued into December 1977 (Fig. 4); however, the temperature 
decrease was almost the same as climatology. At this time, 
the majority of the region had below normal temperatures 
and the anomaly had penetrated 60 m (Fig. 4). From 
Becember 1977 to January 1978 (Fig. 5), the rate of tempera- 


ture decrease (Fig. 6) slowed to half that of climatology 
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in the northeast region, whereas the rate increased to 
1.5 times climatology in the southeast region as the anomaly 
expanded southward. 

The anomaly continued to increase in amplitude and 
Spatial extent between January and March 1978. Temperature 
decreases became more erratic over the area during this 
period. The vertical structure of CA 12 during January 
1978 and March 1978 is shown in Figs. 5 and 7. The cold 
anomaly penetrated to 120 m in January 1978. A warm 
anomaly, extending to 120 m, was found along the eastern 
boundary in January 1978. This resulted in an anomaly 
pattern with a north-south orientation. During February 
(not shown) and March 1978 this warm anomaly extended 
westward in the vicinity of 44°N, 140°W. Thus CA 12 had 
a maximum amplitude in the southeast region. During 
March 1978, CA 12 reached maximum intensity (-1.9l1 °C at 
36°N, 150°W) and areal extent. The western anomaly pattern 
Mmarethe vicinity of 38°N, 170°E in March 1978 was due to the 
formation of CA 18 (Elsberry, 1983) which persisted through- 
out the remainder of the existence of CA 12. 

The vertical structure for April through June 1978 is 
shown in Figs. 8 to 10. The anomaly amplitude diminished 
meeme March to April 1978. During this period, the tempera- 
ture in the vicinity of the central point increased at a 
rate (not shown) of 1.5 times the climatological rate as 


this area underwent spring transition. However, the 
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temperature in the eastern portion of the region remained 
below normal through June 1978 (Figs. 8-10). Through that 
time, the anomaly center was maintained but the pattern 
changed shape. From April to May 1978 the surface features 
of the anomaly became more diffuse, but the anomaly ampli- 
tude at 60 m showed little change. The most rapid surface 
temperature increase took place during these two months 
(Fig. 11); it was associated with the very late spring 
transition throughout the remainder of this region 
(Elsberry, 1983). In the region of CA 12, the temperature 
increase during May to June 1978 was much the same as the 
expected climatological increase (Fig. 12). However, the 
region of CA 18 increased in temperature at twice the 
climatological rate. In May 1978 the surface signature of 
CA 12 was absent, although there was a 10° longitude west- 
ward displacement of minimum temperatures at 60 and 120 m. 
By June 1978 (Fig. 10), the cold anomaly 12 became diffuse 
at 60 and 120 m, but interestingly, it intensified again 
at the surface. In June 1978, there was substantial vertical 
tilt between the surface and 60 m but not between 60 and 
120 m. After July 1978 (not shown), the region was generally 
cold throughout the remainder of the year and periods with 
above normal temperatures did not persist (Elsberry, 1983). 
In summary, CA 12 was a long-duration, large-scale 
event with multiple periods of amplification followed by 


Sustaining periods. The location and minimum temperature 
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for CA 12 each month is given in Table 1. CA 12 formed 
between October and November 1977 at about 44°N, 145°W and 
was only sustained by the expected cooling during the 
southward expansion in December 1977. In January 1978 the 
anomaly progressively expanded southward. The anomaly 
reached its maximum amplitude and areal extent in March 
1978. By April 1978, the anomaly reached its southern 
most location (32°N, 155°W) and had started to diminish in 
Size and amplitude, although portions of the regions 
remained colder than climatology for several months. In 
May 1978, the surface signature of the anomaly disappeared. 
In June 1978, CA 12 was displaced northward to about 38°N, 
155°W and the signature of the anomaly had all but disap- 
peared at depth and had reappeared at the surface. This 
cold anomaly was clearly a major climatolgoical event. 
These changes in the analyzed strength of CA 12 (Table 1) 
could be due to several effects: (1) atmospheric forcing 
pulses due to localized atmospheric effects; (2) oceanic 
response to seasonal climatic changes; (3) anomaly-anomaly 
interactions; and (4) advection effects resulting from 
Ekman and/or Sverdrup mass transport. An effect of the 
above physical mechanisms is that CA 12 could be interpreted 
as being composed of several consecutive cold anomalies 
that have monthly life spans and appear in a region that 
1s anomalously cold. It is not possible to separate 


completely the above effects with the available data. 
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TABLE l 


Comparison of analyzed CA 12 and model MA 12 development for 
each month. The locations and temperatures given are the 
location of the lowest temperature in CA 12 and MA 12 for 
each month. 


ANALYSIS MODEL 

Meme Year Location Temperature Location Temperature 
Nov m977 44°N,145°W -1.63°C 46°N,140°W aoe C 
Dec 1977 40°N,145°W <-1.61°C (1) 46°N,140°W HIE ASS) Oe 

SON So ay =—U48 °C 
Jan ie7s 38°N, 145°W -1.36°C (1) 46°N,140°W See 

(2334-55 W =2.14°C 
Feb 1978 34°N,145°W <-1.76°C 34°N,150°W = 22. 23'2C 
Mar mye 36°N,150°W -1.91°C 34°N,150°W ele 
Apr m72 32°N,155°W -1.85°C SOs es SW =o) LC) e 
May iwc s N, VS0°CW* —1.31°C C0 he ao WwW =2,65 C 
Jun ica sa-N,;ios W —-1.65°C S6-N7 lS 02W -3.79°C 
beet 60 m 
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There may be some tilt present in the vertical structure 
of the anomalies. These tilts will not be examined here 
because the objective analysis of the temperature field 
is done independently at each level, and the horizontal 
resolution of these analyses 1S not appropriate. The 


next section describes the model hindcast of CA 12. 
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Temperature (°C) anomaly CA 12 at surface 
during (A) October 1977 and (B) November 1977. 
The horizontal temperature analysis for each 
month is centered on the 15th of that month. 
Negative (dashed) lines represent regions with 
temperatures less than climatology, zero (heavy 
solid) same as climatology and positive (light 
solid) temperatures greater than climatology. 
The interval is 0.5 °C. Cross hatched areas 
have insufficient data for analysis. Dot marks 
the ultimate location of maximum intensity or 
Ceontadiapormee(s6°N, 250°W) for CA 12. 
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Figure 3. The ratio of change in analyzed temperatures 
(°C) from 15 October to 15 November 1977 to 
the climatological temperature change for 
the same period. Values greater than 1.0 
during this period indicate above-average 
temperature decreases. 
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Temperature (°C) anomaly CA 12 during 


(A) surface and (B) 60m. 


December 1977 at 


Figure 4. 
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Figure 5. Temperature (°C) anomaly CA 12 during January 
ie eeacmaApeesiirtace, (8B) 60 m and (C) 120 m. 
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Figure 6. Similar to Fig. 3 except for the period from 
15 December 1977 to 15 January 1978. 
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Beagure 8. Similar to Fig. 5 except for April 1978. 
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Figure 10. 
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meeore 11. Similar to Fig. 3 except from 15 April to 


15 May 1978. Values greater than 1.0 indicate 
greater temperature changes than the normal 
monthly change. Negative (dashed) values 
indicate temperature decrease during a period 
in which seasonal warming is taking place. 
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Pmeure 12. Similar to Fig. 3 except from 15 May to 


15 June 1978. 
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B. MODEL HINDCAST 

The Garwood model was initialized with the ocean tem- 
perature analysis from 15 October 1977 and integrated 
through 18 June 1978 to cover the life span of CA 12. As 
may be seen in Fig. 13, the model does a rather remarkable 
job in predicting the mixed layer temperature through 
March 1978 during the deepening and maintenance phases of 
the anomaly at the central point (36.0°N, 150.0°W). From 
April through June 1978 however, the model-predicted mixed 
layer temperatures remained too low. 

Recall that the model was initialized Caer the October 
1977 fields which showed the region covered by above average 
temperatures (Fig. 2A). By November 1977 (Fig. 14), the 
mode. developed a cold anomaly on the surface in the north- 
east region. The model anomaly (MA) underwent rapid 
temperature decreases (Fig. 15) from October to November 
1977 as did CA 12. In the northeast region, the analyzed 
surface temperature decrease was greater than 1.5 times the 
Bekiato logical decrease (Fig. 3), while the model surface 
temperature decrease was greater than twice the climatologi- 
cal decrease (Fig. 15). From November to December 1977, 
the model sustained the anomaly by a cooling rate much the 
same aS climatology. By this time the majority of the 
region had below normal temperatures. The model-derived 
anomaly penetrates to 60 m at the central point in December 


1977 (Fig. 16), which is one month earlier than analyzed 
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(Figs. 4 and 5). That is, the model decreases in mixed 
layer temperature (and increases in depth) are slightly 
larger than observed in CA 12 from October to November 1977. 
The model anomaly in the northeast (MA 12N) at 46°N, 140°W 
is in the locality where CA 12 first appeared. MA 12N 
remains in the same location for three months, and disap- 
pears aS a Separate center by February 1978. 

In December 1978, another cold model anomaly (MA 12) 
center appears on the surface in the southeast (-1.48°C at 
30°N, 155°W). Both centers (MA 12 and MA 12N) are con- 
Sidered to be part of the model's representation of CA 12. 
As the southeast center (MA 12) development after December 
1978 closely parallels that of CA 12, the discussion below 
will focus on MA 12. From December 1977 to January 1978 
(Fig. 17), the predicted rate of temperature decrease 
(Fig. 18) slowed to that of climatology in the northeast 
region, but stayed at 1.5 times climatology in the southeast 
region as the anomaly expanded southward. This general 
pattern is consistent with the analyzed changes in the 
eastern half of the domain (Fig. 6). 

The vertical structure for January and March 1978 can 
be seen in Figs. 17 and 19. The model predictions are 
noisy at 60 and 120 m, although the features of CA 12 are 
still very discernable. The model prediction takes on a 
weak north-south orientation similar to the analyses. The 


warm anomaly along the eastern boundary is not predicted in 
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January 1978 except at 120 m, where it is correctly oriented 
north-south but has a larger areal extent. By March 1978 
(Fig. 19), this warm anomaly at 120 m has moved eastward 

to along the eastern boundary of the ADS area, where it 
remains through June 1978. Since this warm anomaly has a 
large areal extent in the eastern region at 120 m, MA 12 
does not penetrate to 120 m in the vicinity of 36°N, 150°W 
until February 1978 (not shown), which is one month later 
than analyzed. 

The model does a fairly good job in hindcasting the 
marked changes in shape of the vertical thermal structure 
for the central point of CA 12. The model mixed layer 
temperatures (MLT's) are in close agreement, albeit on the 
cold side, with those analyzed through March 1978. The 
analyzed and predicted vertical temperature profiles for 
October through December 1977 can be seen in Fig. 20A and 
PeBeand those for January through March 1978 in Fig. 21A 
and 21B. The October 1977 profiles are, of course, identical 
as this is the initial time. The predicted mixed layer 
temperatures (MLT's) for November and December 1977 were 
nearly the same as analyzed, while the remainder of the 
vertical profiles approximated the analyses. One exception 
is the higher thermoclime temperatures analyzed during 
November 1977 only. It is during the period from January 
to March 1978 that the model predictions come the closest 


to hindcasting correctly the vertical temperature profiles 


44 





at the central point down to 200 m. The mixed layer depth 
is also largest during this period. The model MLT's for 
January through March 1978 are less than 0.5°C lower than 
the analyzed MLT and the thermal structure below the 

mixed layer is in fair agreement with the analyses. 

The analyzed spring transition at the central point 
occurred between March and April 1978 using the criteria 
described in Chapter I, Section Cl. Determination of the 
mixed layer depth for both the model and analyzed profiles 
was made using the criteria described in Chapter I, 
Section C2. The analyzed March 1978 vertical temperature 
profile (Fig. 21A) appears to have a mixed layer depth (MLD) 
of about 80 m, although the temperature gradient below 
this level is quite small. By April 1978 (Fig. 22A), the 
analyzed MLD has shallowed to about 40 m and the MLT has 
increased by 1.0°C. Using the same criteria as for the 
analyzed spring transition, the model spring transition 
does not occur until between April and May 1978. As may 
be seen in Fig. 22B, the model MLD in April 1978 is 150 m. 
By May 1978, the model MLD has shallow-d to about 5 m and 
the model MLT has increased by 1.0°C. Using the daily 
maximum model MLD and corresponding MLT from Fig. 13, the 
predicted spring transition would be during the latter 
Melt of April 1978. 

Throughout the spring, the model mixed layer temperature 


remained too low (Fig. 13). In the southwest region, the 
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model predictions continued to decrease in temperature 

from March to April 1978, whereas the region actually had 
undergone spring transition and the cold anomaly was 
diminishing. The model-derived vertical structure for 
April through June 1978 can be Seen in Figs. 23 to 25. It 
is during April 1978 that MA 12 reaches maximum intensity 
f—3eL0°C at 30°N, 155°W) and areal extent. The model 
anomaly center in April 1978 is southwest of the analyzed 
central point. Between April and May 1978, the surface 
temperature in the model began to increase at a rate 

(Fig. 26) close to twice the climatological rate, which is 
consistent with the analyzed changes (Fig. ll). This rapid 
increase in temerature did not persist long enough for the 
model mixed layer temperature to increase to the analyzed 
temperatures. In May 1978 (Fig. 24), the predicted anomaly 
is quite intense rather than having somewhat diffuse surface 
features, as was analyzed. The model prediction does 
Maintain the eastern portion at below normal temperatures 
through June 1978. From May to June 1978, the predicted 
rate of temperature increase (Fig. 27) is less than that 

Of climatology, and is approximately equal to the analyzed 
rate of temperature increase (Fig. 12) in the region of 

CA 12. By June 1978, at the central point, the model- 
derived anomaly is still quite prominent, -1.79°C, at 


60 m, whereas at 120 m the anomaly has diminished to only 
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-0.65°C. The model correctly predicts the intensification 
of the surface anomaly at the central point in June 1978 
(Fig. 25); however, the model anomaly iS twice as intense 
as analyzed. 

In summary, the Garwood model does a very notable job 
in hindcasting CA 12 as it developed during the autumn of 
1977 and reached the maximum amplitude and area extent 
during the winter of 1978. The mixed layer temperature 
and depth, and corresponding vertical temperature profiles 
at the center of the anomaly, are fairly close to those 
analyzed. The model spring transition is late compared to 
the analyzed spring transition, and the model mixed layer 
temperatures remain too low, A summary of the locations 
and minimum temperatures during the life spans of both 
CA 12 and MA 12 is provided in Table 1 (Chapter III, 
Section 3B). It is very satisfying to see numerous sSimilari- 
ties during the life span of the cold anomaly. The model 
hindcast of CA 12 is slightly on the cold side as seen in 
Fig. 13. These results show that CA 12 can be accounted 
for the most part by local atmospheric forcing. 

It iS interesting to look at other grid points around 
the central point of CA 12. One grid point examined is 
located to the northeast of CA 12 near the persistent warm 
anomaly along the eastern boundary of the ADS area. The 


model predictions there are not as good, as the model does 
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not hindcast the warm anomaly along the eastern boundary. 
At 38°N, 135°W (Fig. 28), the model mixed layer temperature 
is also consistently too low throughout the integration. 
However, the monthly trends of temperature change are well 
Simulated. At the northwest edge of the anomaly (40°N, 
160°W), themodel prediction (Fig. 29) is much better than 
at the central point of CA 12. Thus, the temperature 
changes in the central area of the ADS region appear to be 
explained by one-dimensional processes, Furthermore, the 
anomalous thermal structure appears to be generated by 
local atmospheric forcing. The hindcast for 34°N, 160°W 
(Fig. 30) is also quite good. At this location, the model 
captures the trend during the autumn-winter temperature 
decrease and also does a good job in the spring, except for 
April 1978. In April, the model mixed layer temperature is 
Significantly colder than the analyzed monthly mixed layer 


memoerature. 
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Figure 14. Model temperature (°C) anomaly MA 12 at 
surface during November 1977. 
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Figure 15. The ratio of the change in model temperatures 
Geimetrom 25 October £O 15 November 1977 to 
the climatological temperature change for the 
Same period. Values greater than 1.0 during 
this period indicate above-average temperature 
decreases. 
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Figure 16. Model temperature (°C) anomaly MA 12 during 
December 1977 at (A) surface and (B) 60 m. 
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than 1.0 indicate greater temperature changes 
than the normal monthly decrease. 
(dashed) values indicate a temperature 
increase during a period in which seasonal 
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See eR ROR DISCUSSION 

A later (earlier) than normal mixed layer spring 
transition is expected to be related to the development of 
lower (higher) than normal sea-surface temperatures. This 
relationship held rather well for the 19-year sample at 
Ocean Weather Station "P" (50°N, 145°W) (Elsberry and 
Garwood, 1978) and also during 1976 and 1977 in the ADS 
meowon (Budd, 1980). If the timing of spring transition 
is not correct (for example, due to incorrect forcing), the 
model predictions after that time will be in error. 

One may observe these errors in the heat content of the 
model vertical temperature profile. The difference in heat 
content relative to 200 m between analyzed and model 


Peotiles, AQ, iS given by: 


0 
AQ. = oS (Tonal ‘2? Tonal $299?) az 
bce (1) 
0 
. honagtes (T og (2) - Tog (220) dz 


where - is the density, c. is the specific heat and z is 
the depth below the surface. The relative error in the 
heat content relative to 200 m is given by: 


AQ 


relative error = 9 —————————— oO... (2) 
Po Wp J oe : aera e2 
-200 
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Trapezoidal integration is used in the determination of the 
heat content. The heat content iS calculated relative to 
the 200 m temperature with the intent of removing some of 
the change in heat content due to vertical displacements 

of the thermocline and also due to horizontal effects. To 
remove completely the change in heat content due to vertical 
displacements, one should integrate with respect to an 
isotherm in the thermocline. However the analyses have 
insufficient resolution to allow such an integration. 
Another motivating factor for calculations relative to 

200 m was that the calculations to determine the correction 
fields to the surface heat flux fields were also relative 
to this depth (Elsberry et al., 1982). 

Values for the analyzed and model hindcast heat content 
relative to 200 m from November 1977 to June 1978 are given 
in Table 2. Also shown in Table 2 are the differences in 
heat content for each month and the relative error in the 
heat content between analyzed and model profiles. The 
model heat content is derived from the temperature profiles 
averaged over five days centered on the 15th of each month. 
Using averaging periods of 10 and 20 days showed little 
change in the heat content of the averaged model tempera- 
eure profiles relative to 200 m. Thus all of the model 
profiles used in this study are based on the five-day 


averaging period. 
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Bee 2 


Analyzed and model hindcast heat content (x104 cal/em*) 
meeertrvye to 200 m at 36.0°N, 150.0°W, their difference, 
and the relative error in heat content above 200 m between 
analyzed and hindcast temperature profiles. 


Month Year Analyzed Hindcast Difference Relative Error 


Nov 77 6233 6.64 203i -0.049 
Dec 1977 6.33 5.79 0.56 0.086 
Jan 1978 4.09 3.95 0.14 0.036 
Feb 1978 Pi Ben -0.04 =), O13 
Mar 1978 2.38 2.48 -0.10 -0.043 
Apr 1978 ome 2.45 0.79 0.244 
May 1978 3.88 2.89 0.99 0.255 
Jun 1978 4.47 3.43 ae lee 


The relative errors in Table 2 show a seasonal trend. 
During winter, the mixed layer is deep and the relative 
error is small. During summer, the mixed layer is shallow 
and the relative error is increased over five-fold. This 
is not Surprising. When the mixed layer is shallow and 
warm during the summer, there is a greater possibility that 
the model profile below the mixed layer will be in error, 
Since it is essentially unchanged once the mixed layer 
shallows. 

The relative error is a function of (1) analysis error; 
(2) model parameterizations; (3) model forcing; and 
Each of these areas will be 


(4) model physical processes. 


examined in more detail in the next four sections. 
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meanalysis Error 
Thesnumser and distribution of ship XBT reports that 


are included in the objective analysis (White and Bernstein, 
1979) determine the accuracy of the fields. One problem 
appears to be a "fair weather bias" because ships tend to 
avoid bad weather. During the winter, the ship tracks are 
further south, and consequently the analysis in the northern 
portions of the ADS grid is extrapolated from data in the 
southern portions. Another problem is that merchant ships 

do not transit the region at equally spaced time intervals. 
There is no assurance that the reports will be well distri- 
buted in space or in time. A candidate for a varied temporal 
distribution of XBT reports is the month of April 1978. 

Since CA 12 is in the southern portion of the ADS grid, 

where the monthly mean wind fields indicate persistent anti- 
cyclonic flow, it is unlikely that ships would avoid the 
region due to bad weather. As seen in Fig. 13, the model 
mixed layer depth has shallowed to about 25 m at the beginning 
of April 1978, however, it is predicted to deepen again to 
about 70 m in the middle of the month and then shallow to 
approximately 10 m in late April. If ship observations were 
taken at the beginning and end of the month when the mixed 
layer was warm and shallow, this would result in an erroneous 
analysis of MLT. As noted in Fig. 30, a comparison of the 
observed and model MLT's to the southwest of the central 
point also shows that the model MLT was lower than the 
observed MLT for April 1978. 
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2. Model Parameterizations 

iets AGeweLOUm COnStants that are critical to the 
parameterization of the physical processes in the Garwood 
model. These constants, M3, D3, © and y, must be evaluated 
empirically (a more complete deScription of these constants 
and the parameterizations is contained in the Appendix). 
The model parameters m3 and P3 are important in predicting 
the seasonal cycle of sea-surface temperature. The values 
of these constants were empirically determined by Gallacher 
et al. (1983) for the ADS region using data between January 
1976 to December 1978. 

The model constants r (fraction of solar flux 
absorbed in first meter) and y (absorption coefficient) 
determine the predicted vertical distribution of heat in 
the water column, and are functions of the turbidity of the 
water. These parameters are moSt important during the 
Spring and summer when the mixed layer depth is small. To 
determine the model Sensitivity to r and y, the model was 
rerun from 15 February to 18 June 1978 with r decreased 
Beem 0.5 to 0.4 and y increased from 0.001 to 0.002. These 
values will lead to more heat being absorbed into the 
surface layers of the column. However, the new mixed layer 
temperature prediction was not a Significant improvement on 
the model results shown in Fig. 13. This was not surprising 
Since the model constants have already been adjusted to be 


as consistent as possible with the atmospheric forcing 
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(Gallacher et al., 1983). Thus, some other reasons/ 
processes are responsible for the analysis-model differences 
in the spring. 
3. Model Forcing 

It is necessary to first determine if the errors are 
due to an initial value problem associated, in this case, 
with the long integration period. That is, the accumulated 
errors in the model forcing or in model constants could have 
led to incorrect profiles prior to the spring period. To 
test the contribution from this source of error, the model 
was initialized with the 15 February 1978 analysis. The 
daily maximum mixed layer depth and corresponding tempera- 
ture plots when the model was initialized in February 1978 
were very nearly the same as in Fig. 13. The relative 
errors in Table 3 for this experiment are similar to those 
in Table 2, which shows that there is minimal improvement 
by initializing the model in February 1978. To determine 
if the model could produce the correct trend after the 
spring transition had occurred, the model was also 
initialized with the 15 April 1978 analysis. Even in this 
case, the model did not predict the correct trend in the 
mixed layer temperature, and the predicted June 1978 sea- 
surface temperatures continued to be too low. 

Since the errors were not due to an initial value 
problem, the actual model forcing from October 1977 through 


June 1978 was examined. As the model predictions were good 
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through March 1978, it is assumed that the forcing through 
that time were reasonably correct. Earlier studies 
(Elsberry et al., 1979; Budd, 1980; Steinder, 1981) found 
that the monthly integrated FNOC surface heating did not 
have a sufficiently large seasonal amplitude. There was 
also a persistent bias toward excessive heat loss to the 
atmosphere especially along the southern boundary of the 
domain. Elsberry et al. (1982) calculated bimonthly heat 
flux correction fields to be added to the FNOC surface 
heat flux fields to reduce this excessive upward flux. 
There was a downward heat flux correction of about 

4 meen /h applied during April-May (Elsberry et al., 
1982). The fraction of the correction fields that should 
be attributed to the solar heat flux was not determined by 
Elsberry et al. (1982), as previous studies (Elsberry et 
al., 1979: Budd, 1980; Steiner, 1981) indicated that the 
solar heat flux fields were reasonably correct. The 
resulting heat fluxes were used in the model forcing for 
feos study. 

As mentioned earlier, the model predicted inadequate 
ocean warming in the spring. This suggests that the most 
likely cause of the observed errors 1S inaccurate atmo- 
Semeric forcing. It is important to point out that the 
predictions during the spring are very sensitive to the 
magnitude of the downward surface heat flux. There is a 


delicate balance between the heat flux tending to form a 
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stable and shallow layer near the surface and the wind 
mixing tending to erode the stable layer to re-establish a 
mixed layer characteristic of winter. Once the shallow 
mixed layer is formed, the greater the downward flux, the 
more rapid the warming, and the less likely that wind 
mixing will erode the layer. Thus, the magnitude of the 
heating has a double effect in the spring transition period. 
As discussed earlier, this was the period when the relative 
errors in heat content relative to 200 m began increasing 
dramatically. 

Budd (1980) found that the time changes in fluxes 
of solar radiation and total heat flux were apparently of 
much less importance than variations in wind speed during 
the change of the mixed layer from a winter to a summer 
regime. Thus the effect of errors, even small ones, in 
wind speed forcing are enhanced during the period of spring 
transition. It is not possible to determine corrections to 
the FNOC wind fields as no indpendent wind reports are 
available. However, one would expect the FNOC analyses to 
be biased toward lower wind speeds, due to the fair weather 
bias. Since the model mixed layer depth is too deep, the 
FNOC winds would have to be too high, which is unlikely. 

To illustrate the sensitivity of the model predictions 
during this period, an additional downward heat flux was 
added to the surface heat flux correction determined by 


Elsberry et al. (1982). Several model runs with additional 
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downward heat flux were made. These runs were initialized 
in February 1978 to reduce computer time. It was found 
that to predict correctly the trend from March to June 1978 
(Fig. 31), an additional downward heat flux of about 

5 PeVem-/h had to be added each time step from 15 March 
to 18 June 1978. With the additional heat, the relative 
errors in heat content dropped dramatically (Table 3) and 
the spring transition occurred during the same period as 
the analysis. The resulting horizontal structure features 
for March and April 1978 (Figs. 32 and 33) were closer to 
the analyzed structure (Figs. 7 and 8) of CA 12 than the 
Original model predictions (Figs. 19 and 23). The addi- 
tional heat also created a small warm anomaly along the 
eastern boundary of the ADS region in March 1978. By 
April 1978, this warm anomaly covered a larger section of 
the northeast than was analyzed, but the anomaly features 
are still closer than the original predictions to the 
analyzed features. The improvement of the model features, 
the presence of a warm model anomaly along the eastern 
boundary and the reduction in relative errors suggests 
that the most likely cause of the errors at the central 
point for March and April 1978 is inaccurate atmospheric 
forcing. However, the downward surface heat flux of 

5 ey om /h did not improve the horizontal temperature 
anomaly structure for May and June 1978 (not shown). Those 


fields had much too strong thermal anomalies which did not 
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resemble the analyzed horizontal structure, although the 
relative errors in heat content were decreased substantially 
(Table 3). This suggests that a correction of about 

5 mcm /h of heat should be added, but not entirely as 
surface heat flux. It is also possible that some other 
reasons or processes are responsible for the analysis-model 


differences found in May and June 1978. 
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TABLE 3 


Relative error in model hindcast heat content relative to 
200 m at the central point 36.0°N, 150.0°W when the model 
was initialized on 15 February 1978 (column labeled 
Hindcast) and when an additional 5 cal/cm2/h was added to 
the downward surface heat flux from 15 March to 18 June 
(last column) 


Hindcast 
Month Meat Hindcast + 5 cal/cm2/h 
Mar 7 =O) een) ape 0-026 
Apr 1978 0255 O233 
May 1978 0.260 0.062 
Jun oS G2236 -0.033 
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4. Model Physical Processes 

This one~dimensional version of the Garwood model 
does not treat advective effects. Of the physical processes 
not included in the model, horizontal advection is the most 
likely to contribute to the erroneous prediction during the 
spring of 1978. To examine this effect it is necessary to 
describe the synoptic situation during the spring. The 
monthly mean streamfunction fields, calculated from the 
FNOC winds over the ADS grid, were used to examine the 
interannual variability in the monthly mean surface pressure 
fields. The streamfunction is generally proportional to and 
parallel to the isobars. | 

During March of 1976, 1977 and 1978 a large anti- 
cyclone was centered near 35°N, 135°W and covered most of 
the eastern North Pacific south of 45°N. This general 
anticyclonic pattern continued through April with the cell 
weakening in both 1976 and 1977, but intensifying in 1978. 
In the northern portions of the ADS area during March and 
April the winds were generally westerly. During May, the 
anticyclone was replaced by a cyclonic pattern, except in 
1978. The anticyclone was anomalously weak over the 
southeast of the region (centered approximately 34°N, 
145°W) during May 1978. In June of all three years, the 
eastern Pacific subtropical anticyclone was re-established 
over the eastern half of the ADS grid. 

Two arguments can be given to explain why the model 
temperatures were too low in May and June 1978. The first 
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argument concentrates on the atmospheric effects of the 
anomalous anticyclonic circulation over the eastern part 
@eecne ADS area during May 1978. With an anticyclonic 
circulation, there would be an increased downward solar 
heat flux due to the clear skies and light winds. Neither 
of these conditions is conducive to turbulent mixing in 
the ocean, so that an increase in downward solar heat flux 
may account for the errors in the spring. 

The second argument examines the anomalous surface 
Ekman and Sverdrup mass transport. The anticyclonic winds 
during May 1978 are southeasterly over the region of the 
central point. This is expected to produce a surface Ekman 
mass transport of warmer water towards the northeast (Since 
the depth-integrated Ekman mass transport is 90° to the 
right of the wind stress) into CA 12. To balance this 
convergence and subsequent downwelling caused by Ekman 
pumping under the atmospheric anticyclone, there must be a 
geostrophic mass transport such that the total Sverdrup 
mass transport balances the negative curl of the wind stress. 
The balancing mass transport is southward bringing colder 
water into the region of CA 12. There is not sufficient 
data to say which mass tranSport is dominant. The observed 
vertical temperature profiles below the mixed layer for 
April, May and June 1978 (Fig. 22A) do show that there was 
warm advection in April, cold advection in May followed 


again by warm advection in June. 
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In summary, the synoptic scale flow over the region 
of CA 12 during April and May 1978 was more anticyclonic 
than during 1976 or 1977. During March and June the synoptic 
scale flow over the region of CA 12 was anticyclonic all 
three years. Either an additional downward heat flux or an 
Ekman transport of warm water would be consistent with the 
synoptic conditions. It is not possible to determine which 


effect 1s more important due to lack of data. 
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fe ue RGE=oCGALE RAPID TRANSITION CASE 


A. COLD ANOMALY 1 DESCRIPTION 

The cold anomaly described below was designated and 
described by Elsberry (1983) as CA 1. CA 1 is a shallow, 
large-scale, rapid transition event that persisted for 
approximately seven months. It was the most intense cold 
anomaly detected in the ADS region during 1976-1979 
(Elsberry, 1983). 

This cold anomaly is especially dramatic as it occurred 
after at least six months (the first available map is 
January 1976) of above-normal temperatures at all depths 
(Elsberry, 1983). The conditions on the surface during 
June 1976 (Fig. 34) show that the majority of the ADS region 
had above-normal temperatures due to a large-scale, warm 
anomaly centered 36°N, 160°W. For orientation purposes, the 
ultimate location of maximum intensity or central point 
(40°N, 165°E) for CA 1 is marked as a dot on the figures. 
The rate of temperature increase (Fig. 35) from June to 
July 1976 (Fig. 36) along the western boundary was less 
than 0.5 times the expected climatological rate. This gave 
Bese tO CA 1, which during July 1976 (Fig. 36) included two 
centers at 42°N, 165°E and 36°N, 165°E with temperatures of 
-2.47PC and -2.18°C respectively. The southern center 
appeared to be relatively deep as it penetrated to at least 


120 m, but this center did not persist. The eastern anomaly 


82 





pattern in the vicinity of 42°N, 140°W in July 1976 (Fig. 36) 
was designated CA 2 by Elsberry (1983). This anomaly per- 
Sisted throughout most of the life of CA 1. The very slow 
rate of temperature increase relative to climatology 
jeg. 3/7) continued into August 1976 (Fig. 38). It was 
during August 1976 that CA 1 reached maximum intensity 
(-4.41°C at 40°N, 165°E) and areal extent. At this time 
CA 1 eradicated the surface Signature of a warm anomaly 
which was located near 42°N, 175°E during July 1976. Even 
though this cold anomaly is very intense, it is interesting 
that it was a very shallow feature. As may be seen in 
Fig. 38, CA 2 also reached maximum intensity (-1.76°C at 
42°N, 140°W) and areal extent during August 1976. 

Summer warming continued erratically from August to 
September 1976 (Fig. 39) over the region of CA Wewneress 
a cooling trend would have been expected based on clima- 
tology. By September 1976 (Fig. 39), CA 1 began to diminish. 
The northern section of CA 1 was sustained, although the 
Senter of -2.32°C was displaced to 42°N, 170°E. The 
southern section diminished during September 1976, perhaps 
in response to the arrival of an intense warm anomaly at 
60 m and below. 

The expected fall cooling occurred between September 
and October 1976 (Fig. 40) over the entire grid as the 
rates of temperature decrease (Fig. 41) were of the same 


Seder as climatology. During October 1976 (Fig. 40),-CA l 
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took on a two-cell appearance with centers at 42°N, 180° 
and 38°N, 165°E and with temperature anomalies of -1.81°C 
and -1.35°C respectively. The western center extended 
below 120 m. Between these two centers was a strong warm 
anomaly at 60 m. From October to November 1976 (Fig. 42), 
the rate of temperature decrease (Fig. 43) over the entire 
grid was again much the same as climatology. Thus the 
anomaly is simply being maintained during September to 
November rather than being Significantly strengthened or 
diminished. In November 1976 (Fig. 42), CA 1 increased in 
intensity with depth and extended through 120 m. By 
December 1976, CA 1 was greatly diminished and had a two- 
_cell appearance through 60 m, with the strongest signature 
at 120 m (Fig. 44). Above-normal temperatures were present 
over much of the ADS region in December 1976 (Fig. 44), as 
the rates of temperature decrease (not shown) were slightly 
less than climatology. The cold anomaly was not discernible 
after December 1976. 

CA 1 often had a two-cell appearance (Table 4). These 
changes in the analyzed strength of CA 1 may be the result 
of the physical mechanisms described in Chapter 3, Section A. 
In November 1976 at 36°N, 165°E and in December 1976 at 
36°N, 160°E, the strength of the anomaly increases with 
depth. There are several explanations for the decrease in 
anomalous temperature with depth. A possible explanation 


is that there was surface Ekman mass transport of warmer 
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water into the region. Based on the monthly mean stream- 
function fields for 1976, 1977 and 1978, there was anomalous 
westerly flow over the region of CA 1 in November and 
December 1976. This flow would not result in the transport 
of warmer surface water into the region by surface Ekman 
mass transport. Instead, it would tend to increase the 
surface intensity of CA 1. Another possibility is cold 
advection at depth. However, a complete description of 

the cause of this anomalous decrease in temperature with 
depth is not possible due to lack of data and grid 
resolution, 

In summary, CA 1 was a shallow, long-duration, rapid 
transition event. CA 1 formed between June and July 1976 
in the west central section of the ADS domain with two 
centers of low temperature (Table 4). Anomaly development 
waS rapid over a region in which the temperature increased 
at a much slower rate than expected from climatology. The 
cold anomaly reached maximum amplitude and areal extent 
during August 1976. The intensity of this anomaly was 
truly remarkable considering that it occurred during the 
summer. In September 1976, the anomaly started to diminish 
in size and amplitude. Up through September 1976, the 
center of lowest temperature tended to be in the vicinity 
of 41°N, 163°E. As CA 1 diminished, the centers of lowest 
temperature erratically moved southward. The anomaly 


persisted above 60 m until October 1976. Portions of the 
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regions remained colder than climatology through December 
1976. The next section describes the model hindcast of 


eA l. 
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TABLE 4 


CA 1 development for each month during 1976. The location(s) 
and temperature(s) are for the lowest temperatures in region 
@: CA 1 for that month. 


Month Location Temperature 
Jul (aee2oN, 165° = 2247 °C 
(2S o oN 165 7 E =2..Le ue 
Aug AON LoS: & -4.41°C 
Sept a2 aN nO & =2e52ae 
Oct Cees ecN, 165° =1755 .¢ 
(2) 42°N, 180° aap ae 
Nov BON Loo 7 = On 2 +e 
Dec Mpc aN 165-5 =o. .C 
CZ) o4°N 7 160° E -0.76°C 
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Figure 34. 


migure 35. 
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L60E 170E 180 170W 160W 150W 140W 130W 
Temperature (°C) anomaly CA 1 at surface during 
June 1976. The horizontal temperature analysis 
for each month is centered on the 15th of that 
month. Negative (dashed) lines represent 
regions with temperatures less than climatology, 
zero (heavy solid) same as climatology and 
poSitive (light solid) temperatures greater than 
Glamacology. The interval is 0.5°C. Cross 
hatched areas have insufficient data for 
analysis. Dot marks the ultimate location of 
maximum intensity of central point (40°N, 

iG) for CA 1. 
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The ratio of change in analyzed temperatures 
(°C) from 15 June to 15 July 1976 to the same 
climatological temperature change for the 

Same period. Values greater than 1.0 during 
this period indicate above-average temperature 
increases. 
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Similar to Fig. 35 except for period from 
Peowoweeco 15 August 1976. Values greater 
than 1.0 indicate greater temperature changes 
than the normal monthly seasonal increase. 


90 





\2 
O 


OOOO; 


<% 


5 


we 


24 
een 


= 
oO 
co 


THI KN 





40N 





SON 


170W 160W 15OW 140W 130K 


180 


[60E 170€ 


36 except for August 1976 


Similar to Fig. 


Preure 38. 


91 








Cy 
LW 
OR 
scat 


: 


RY oe 
Sc 
ef 
2 


AAAS 


M 
x5 
NITRO 


. 


40N 





SON 


L70W 160W 150W 140W {130W 


180 


Reine! / OE 


Similar to Fig. 


Figure 39. 


36 except for September 1976. 


Bg 





8 


\ 


o 





ie 


170W 160W 1SOW 140W 130 


Similar to Fig. 36 except for October 1976. 


180 


Roiew ly Vic 


Figure 40. 


7 





SON 


45N 


40N 


Figure 41. 





OW 


Similar to Fig. 35 except for the period 

15 September to 15 October 1976. Values 
greater than 1.0 indicate greater temperature 
decreases than the normal monthly change. 
Negative (dashed) values indicate a 
temperature increase during a period in which 
seasonal cooling is taking place; except, in 
the lower right hand corner where the dashed 
lines indicate a temperature decrease in an 
area where climatology iS warming. 
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Figure 43. Similar to Fig. 35 except for the period 
15 October to 15 November 1976. Values 
greater than 1.0 indicate greater temperature 
decreases than the normal monthly seasonal 
change. 
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B. MODEL HINDCAST 

The Garwood model was initialized with the TRANSPAC 
temperature analysis from 15 June 1976 and integrated 
through 18 December 1976 to cover the life span of CA l. 
Instead of hindcasting the extremely intense CAl, the first 
model prediction created an extremely shallow and warm 
anomaly. The model results for August 1976 at the surface 
can be seen in Fig. 45. The daily maximum model mixed layer 
depth and corresponding mixed layer temperatures (not shown) 
at the central point (40°N, 165°E) of CA 1 shows that the 
temperatures are over 15°C higher than analyzed. By 
December 1976, the mixed layer temperature is about 5°C 
higher than analyzed. As shown in Fig. 45, these high 
surface temperatures are part of a zonal surface temperature 
pattern with temperatures much too high in the regions south 
of 42°N and too low in the northern regions. The predicted 
surface temperatures are particularly high to the west of 
180° longitude. The next section will investigate the 


reasons and processes that could produce erroneous results. 
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Figure 45. Model temperature (°C) anomaly MA 1 at the 
surface during August 1976. The model used 
the surface heat flux correction relative to 
200 m as determined by Elsberry et al. (1982). 
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See eR ROR DISCUSSION 

The model anomaly field bore little resemblance to that 
analyzed, and did not produce a cold anomaly in the region 
of CA l. The poor model results could be a function of 
(1) data quality; (2) model physics; and (3) model forcing. 
Each of these areas will be examined in more detail in the 
next three sections. 

1. Data Quality 

The original objective analysis procedure used by 
White and Bernstein (1979) generates values at all grid 
points regardless of whether any observations are present 
in the vicinity of the grid points. In a data-sparse region, 
a faulty XBT could significantly degrade the quality of the 
monthly objectively analyzed values at the surrounding grid 
points and cause a temperature bias with depth. However, 
one would expect that the objective analysis technique 
would tend to eliminate such a bias if other correct profiles 
are in that region. 

CA 1 is the most intense anomaly observed by Elsberry 
meees) during 1976-1979 (Fig. 46). CA 1 is about 2.5°C 
colder than the next coldest anomaly. It is possible that 
the extreme intensity of CA 1 is due to erroneous data. 
However, neither the XBT profiles nor the data distribution 
Mmeps for the life span of CA 1 are available to determine 


if this factor can explain part of the model discrepancy. 
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Histogram showing the occurrence of cold 
anomalies during the period January 1976 to 
December 1979. CA 1 is on the extreme right 
in the histogram (from Elsberry, 1983). 


tOee 





2. Model Physics 

Because Garwood (19/7) found that shear production 
Significantly affected entrainment mixing in relatively 
limited conditions, the version of the model used here does 
not contain equations for mean momentum or a term for shear 
production in the turbulent kinetic energy equation. 
However, there iS some evidence as reported by Price (1981), 
Kraus (1981) and Martin (1982) for the importance of shear 
aS a mixing mechanism. Martin (1983) suggests that the 
Garwood model would benefit from the inclusion of shear 
production due to the mean current, in that augmented shear 
production should prevent overly shallow mixed layers from 
persisting. In the Garwood model, the time derivative 
terms in the entrainment equation are dropped when the layer 
is shallowing, and the new mixed layer depth is solved for 
algebraically. The absence of entrainment shear production 
may account for some of the model tendency toward too high 
Sea-surface temperatures during the spring and summer. 

The parameterization of the absorption of solar 
radiation in the vertical may be another source of error. 
A discussion of the parameterization of the flux of solar 
radiation, Qa, in the Garwood model is provided in the 
appendix. Fifty percent of the incoming solar radiation is 
absorbed in the upper meter of the ocean. The fraction 
that penetrates to greater depths can have a significant 


effect on the development of the upper ocean thermal 
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structure. In spring and summer, solar radiation provides 

a direct means of heating the seasonal thermocline below 

the mixed layer. During cloudy periods and storms, the 
reduction of the solar radiation means less heat will be 
directly absorbed in the mixed layer, which reduces the 
thermal gradient below the first meter, and allows for 
deeper mixing. The turbidity of the water also affects the 
distribution of solar radiation in the vertical. An increase 
in turbidity (not included in the model) Curing the spring 
and summer could increase the sea-surface temperature and 
the stratification below the mixed layer, while decreasing 
the mixed layer depth and the warming of pe ocean below the 
mixed layer. 

CA 1 might also be caused by the -.advection of a cold 
ocean eddy into the region. White and Bernstein (1979) 
emphasize that the baroclinic eddy activity west of 175°W 
was considerably greater than to the east, where the large- 
scale variability dominates. The transition near the 
Emperor seamount chain between the strong and the weak eddy 
regions was quite abrupt, of the order of 5-10° of longitude 
(Bernstein and White, 1977). Eddies south of the Kuroshio 
extension would be cold, which agrees with the anomalous 
temperatures in the region of CA l. The decay of a cold 
eddy is also consistent with the analyzed temperature 
Structure. That is, a cold core eddy could have arrived 


in the region of CA 1 between July and August 1976 and as 
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the eddy decays the region would be increasing in temperature 
@@emstill be colder than climatology. The transitional 

speed of CA 1 can not be determined accurately using the 
monthly position(s) of maximum intensity (Table 4). The 

Mera Size of 2° latitude by 5° longitude is too coarse to 
describe the evolution of mesoscale events. Using the -2.0°C 
isotherm in Fig. 37 to define the edge of CA 1 in August 

1976 gives a radius of about 650 km. This radius is approxi- 
mately 13 times the baroclinic Rossby radius of deformation. 
Thus if CA 1 was a cold eddy, one would expect a large 
vertical extent. However, this is not seen as CA 1] in 

August 1976 (Fig. 38) is completely above 60 m. Since the 
vertical extent is so small and the horizontal scale so large, 


e 


it is unlikely that CA 1 is one large eddy. 
pee model Forcing 

It is necessary to first determine if the errors are 
due to an initial value problem associated with initializing 
after the spring transition when the mixed layer is shallow 
and warm. According to Martin (1983), errors due to model 
initialization are minimized by initializing during the 
winter when the mixed layer is deep. After the mixed layer 
shallows in the spring and the seasonal thermocline begins 
to form, the model results are fairly independent of the 
meetial conditions. To test the contribution from this 
source of error, the model wasS initialized with the 


15 February 1976 analysis and integrated through 
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membPecember 1976. The results of this test showed that the 
errors are not due to an initial value problem as again the 
mixed layer was too shallow and too warm throughout the 
integration period, and especially during the summer months. 
Since the first model predictions with the bimonthly 
heat flux correction fields determined by Elsberry et al. 
(1982) were too warm in the region of CA 1, the model was 
run without the surface heat flux corrections. The inte- 
gration period was from 15 June to 18 December 1976. The 
daily maximum model mixed layer depths and corresponding 
mixed layer temperatures at the central point (Fig. 47) 
were substantially improved over the first prediction. The 
predicted mixed layer temperature anomaly for August 1976 
was within 1°C of the -4.41°C value. However, large negative 
temperature errors occurred in the remainder of the ADS 
domain. This can be clearly seen in the model hindcast for 
the surface during August 1976 (Fig. 48). The errors were 
largest along the southern boundary of the ADS region. 
Plots (not shown) of daily maximum mixed layer depth and 
corresponding temperature at 38°N, 160°W and at 42°N, 135°W 
showed that a seasonal thermocline was not established. 
These results agree with earlier studies (Elsberry et al., 
m979; Budd, 1980; Steiner, 1981) which found that the 
monthly surface heating from FNOC does not have a suffi- 
ciently large seasonal amplitude. These studies also found 


a persistent bias toward excessive heat loss along the 
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southern boundary of the domain. Thus a downward surface 
heat flux correction is definitely needed in the region, 
Sven if it is not at the location of CA 1. Consequently, 
it was decided to examine carefully the June-August surface 
heat flux correction fields as determined by Elsberry et 
mee (1982). 

For comparison, the model results with the correction 
field during the summer of 1978 were examined. The model 
was initialized with the 15 June 1978 analysis and inte- 
grated through 18 August 1978 to illustrate the effect of 
the June-August surface heat flux correction fields. If 
the model did not predict the anomalous conditions by 
August, then it would not get the anomalous conditions after 
August correct either. Model runs were made with and without 
surface heat flux correction fields. From the model inte- 
gration without the correction fields, it was very apparent 
that as for 1976, a downward heat flux correction is needed 
to correct a bias of excessive upward heat flux over the 
ADS region. The model integrations with the correction 
fields predicted mixed layer temperatures that were too 
cold in 1978 compared to being too warm in 1976. fThus, 
the model integrations suggest that the corrected heat 
flux underestimates the heat flux required for June to 
August 1978, and overestimates that required for 1976 in 


the region of CA l. 
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The analyzed vertical temperature profiles at 
40°N, 165°E for June, July and August of 1976, 1977 and 
1978 were examined. The evolution of the profiles during 
two (1977 and 1978) of the three years were very similar, 
while the period of CA 1 development during 1976 was quite 
different. The analyzed vertical temperature profiles for 
August of 1977 and 1978 were characterized by shallow, 
warm (<5 m, about 24°C) mixed layers. Even though the 
mixed layer was also shallow in August 1976 (<5 m), the 
temperature was about 7°C colder. 

A comparison of the model forcing without the 
surface heat flux correction from June to December for the 
three years snowed that 1976 had noticeably less (about 
10-20 Bat/em-/h) total upward flux than the other two years 
over most of the period. This is consistent en the 
analyzed mixed layer temperature for August 1976 being much 
lower than during the other two years. There were no 
dramatic differences in the wind speed or solar radiation 
fields among the three June-August periods. 

In summary, the surface heat flux forcing is 
Getinitely a potential error source. Although the heat 
correction field derived by Elsberry et al. (1982) may be 
appropriate to the normal large-scale evolution, it 1s not 
appropriate to model predictions of the development of 
CA 1 from June to August 1976. Thus a new surface heat flux 
correction field for the period June-August 1976 is 
described in the next section. 
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Figure 48. 
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Pore viseED HEAT FLUX CORRECTION 

CA 1 occurs during a period of ocean warming (4H > 0) 
and net downward ( but less than expected from climatology) 
surface heat flux (Qn > 0) which produces excessive warming 
{AH - Qn > 0). Elsberry et al. (1982) determined that this 
type of error would occur 26.8% of the time using the 
corrected total heat flux. The heat flux correction fields 
for the June-to-August period as determined by Elsberry et 
al. (1982) are based on only three years. Their basic 
assumption 1s that the error in the surface heat flux is 
systematic and that an average over the three years will 
give a stable estimate of the required correction. However, 
it was shown above that the changes in heat content during 
1976 were markedly different. Therefore a new heat flux 
correction field is derived specifically for the 15 June to 
18 August 1976 period. The purpose is to determine if the 
model can predict the correct vertical distribution of the 
heat if an improved estimate of the total surface heat flux 
1s provided. 

The change in oceanic heat content relative to 200 m 
1s 


fH ee Hites At) = Hey (3) 


where At = two months (June - August 1976) and H relative 


to 200 m is given by 


0 
FC) Gee 20.0 laa 5 (9) 
200 
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It iS expected that AH will be positive during this period 
of net downward heat flux (roughly between April and 
September ). A two-month interval was used to be consistent 
with the heat correction field derived by Elsberry et al. 
@eea2z)- They found that bimonthly correction fields were 
Simple to apply and representative of the seasonal variation. 
This two-month interval is also the time during which CA 1 
reached maximum intensity and spatial extent. 

The following budget equation is assumed: 


At 
AH = it Qdt + Residual = Q At + Residual . (5) 


0 Jk 
The integral sign indicates that the air-sea flux is summed 
over the same time interval that AH is evaluated. Qn 1s 
defined as the total surface heat flux and can be expressed 
as the time integral of: 
Where the subscripts s, b, h and 1 refer to solar, back, 
sensible and latent heat fluxes through the sea surface. 
QsF is surface heat flux. The residual term includes non- 
local physical effects (especially horizontal advection) and 
the errors in estimating the heat content changes and the 


Surface fluxes. This budget equation assumes that vertical 


processes dominate the horizontal, and for the space and 
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time scales in this study, that the local change in heat 
content over the given period should be balanced by the 
vertical flux of heat at the air-sea interface. 

The analyzed temperature fields and the model physics 
were assumed correct. Using predicted temperature fields 
for August 1976 determined without any heat corrections and 
the analyzed temperature fields for August 1976, a new 
bimonthly surface heat flux correction relative to 200 m 


at each grid point is 


correction = Se moctnua) - AH nal (Aug) ?/4* 


(7) 


( )fat 


Hod (Aug) . anal (Aug) 


The correction field (Fig. 49) with units of Seu Pe 
was added to the FNOC surface heat flux each hour in the 
model integration. The solar radiation fields are assumed 
to be correct. Negative values in the resulting correction 
field (Fig. 49) indicate that the upward heat flux is to be 
reduced by the amount shown. It can be seen that Hn and 


od 
H are approximately equal at the central point (40°N, 


anal 
165°E), so that only a small correction is required. 

The model was initialized with the 15 June 1976 analysis 
and integrated through 18 August 1976 (Fig. 50). The 
intensity of CA 1 is correctly hindcast, although the 
gradients around the anomaly are too large. There is also 


a warm anomaly along 170°E which is not present in the 


analyzed anomaly field (Fig. 38A). The shape of the 
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anomalous warm region in the mid-Southern section of the 
ADS domain iS in agreement with that analyzed, albeit with 
much too high temperatures. It may also be seen in Fig. 50 
that there is some weakening of the anomalous features 
through the middle of the ADS domain. The model does 
hindcast the eastern region as cold, which agrees with the 
analysis. 

At the central point of CA 1 (Fig. 51), there is very 
good agreement between the monthly objective analyzed mixed 
layer temperature and the temperature corresponding to the 
daily maximum mixed layer depth. The good agreement at this 
point was expected from Fig. 47 since the new correction is 
almost zero. Major improvements due to the new correction 
field relative to either the original (Fig. 45) or no 
(Fig. 48) correction are found through most of the domain. 

To determine if further improvements in the surface 
anomaly field could be made, two other correction fields 
were derived. As the summer mixed layer is shallow and 
warm, the local atmospheric effects are only felt in the 
upper 100 m. Thus a correction field based on the heat 
content changes relative to 100 m was produced (Fig. 52). 
The model was then initialized with the 15 June analysis 
and integrated with this correction field added to the 
FNOC surface heat flux fields through 18 August 1976. The 
resulting surface anomaly field (Fig. 53) is less representa- 


tive of that analyzed than in Fig. 50. The eastern, middle 


bis 





and southern sections were not hindcast well. In the 
region of CA 1, the temperature of the predicted anomaly 
is too low and the anomaly center is displaced northward. 

Since the anomaly in August 1976 is present above 60 Mm, 
it is possible that the heat flux correction field should 
be based on the heat content changes in the upper 50 m. 
Therefore, a correction field relative to 50 m (Fig. 54) was 
also produced. The model was then integrated over the same 
time period as previously. The resulting surface anomaly 
field (Fig. 55) is a much poorer representation of the 
analyzed field (Fig. 38A) than the predicted fields with 
surface heat flux correction fields relative to 100 and 200 m. 
One region of improvement in Fig. 55 is in the eastern section, 
although the improvement is small relative to Fig. 50. 

The pattern correlation was calculated between the 
analyzed and predicted temperature changes between 15 June 
and 15 August 1976. These calculations are over an area 
bounded by 34°N, 44°N, 170°E and 135°W. The pattern correla- 
tion for the model integration with the surface heat flux 
correction relative to 200 m was 0.49, while the pattern 
correlation for the model integration with the correction 
field relative to 100 m was -0.05. Based on these and the 
previous results, the new correction field relative to 200 
mwas judged to be better. The improved results with 
the larger depth suggest that this calculation is more 
effective in removing some of the effects of horizontal 


advection and changes in heat content due to vertical 
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displacements of the thermocline. Making the correc- 
tions relative to 50 and 100 m may intertwine these 
effects with the effects of local atmospheric forcing. 

In summary, if the best possible estimate of the surface 
heat flux 1S provided, the model predictions were improved, 
although they were only a fair representation of the 
analyzed August 1976 fields. One possible source of error 
is the assumption that the total correction should be 
applied to the surface heat flux. Perhaps a fraction of 
the correction field should be applied to the solar heat 
flux. The areas of erroneous high temperatures in the 
model predictions may also be due to the absence of shear 
production due to the mean current. It is possible that 
non-local atmospheric forcing and events may contribute to 
the poor model performance. Given the available data, it 
iS not possible to separate the effects of three-dimensional 
processes, poor forcing, inaccurate ocean analyses and 
incomplete model physics and parameterizations. 

The very strong anomalous conditions associated with 
CA 1 are at least partly due to the anomalous atmospheric 
forcing during the period. Unfortunately, the magnitude of 
the anomalous forcing appears to be of the same order of 
magnitude as the bias in the FNOC surface heat flux fields. 
When the uncertainties in the heat correction field are 
large, the ability of the model to predict anomalous ocean 


conditions cannot be satisfactorily tested. These tests 
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with a single case would suggest that improvements in the 
FNOC surface heat flux are required. It is hoped that the 
new Navy Operational Global Atmospheric Prediction System 


will provide these improved fields. 
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Correction field (cal/cm*/h) relative to 

BUG nimeeG 15 June to 18 August 1976 to be 
applied to the FNOC surface heat flux fields. 
Negative values indicate that the upward heat 
flux is to be reduced by the amount shown. 
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Model temperature (°C) anomaly at the surface 
during August 1976 resulting from the 
correction field relative to 200 m being 
applied to the FNOC surface heat flux fields. 
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Figure 51. Daily maximum model mixed layer depth and 
corresponding model mixed layer temperature 
(MLT) resulting from the new correction field 
relative to 200 m being applied, compared _to 
objectively analyzed MLT (denoted by x on the 
15th of each month) at the center of CA l 
AVON, 65 By from 12° June to 18 August 1976. 
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meggre 52. Similar to Fig. 49 except for correction 
field relative to 100 m. 
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Figure 53. Similar to Fig. 50 except for the application 
of the correction field relative to 100 m. 
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meee 54. Similar to Fig. 49 except for correction 
field relative to 50 m. 
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of the correction field relative to 50m. 
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V. CONCLUSIONS 


The Garwood model does a very commendable job of hind- 
casting the large-scale long duration anomaly CA 12 through 
Emewautumn Of 1977 and winter of 1978. Errors in the model 
prediction increased substantially in the spring, following 
the analyzed spring transition between 15 March and 15 April 
1978. It was found that an additional downward surface heat 
mint Of 5 ee ema during this period produced improved 
predictions. With the additional heat flux correction the 
model spring transition occurred during the same period as 
the analysis. For March and April 1978, the improvements of 
model features, the presence of a warm anomaly along the 
Saeeern ADS region boundary and the reduction in relative 
heat content errors, suggest that the most likely cause of 
errors at the central point for these two months is 
inaccurate atmospheric forcing. However, the additional 
downward surface heat flux did not improve the horizontal 
temperature anomaly structure for May and June 1978. Those 
fields had much too strong thermal anomalies which did not 
resemble the analyzed horizontal structure, although the 
relative errors in heat content were decreased substantially. 
This suggests that a correction of about 5 ow Bae hs of heat 
should be added, but not entirely as a surface heat flux. 


In May and June 1978, there was anomalous anticyclonic flow 
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Over the region of CA l. Thus, either an additional down- 
ward solar heat flux and/or Ekman transport of warm water 
could be the process by which the additional heat needs to 
be distributed in the model. 

These results verify the hypothesis that CA 12 was 
primarily generated by the vertical mixing processes during 
the autumn of 1977 and the winter of 1978. That is, the 
mixed layer temperature change is related to one-dimensional 
processes to the first order. The usefulness of the heat 
flux correction fields, derived by Elsberry et al. (1982), 
for FNOC surface heat flux estimates in North Bacaet Ocean 
predictions was demonstrated during the autumn of 1977 and 
the winter of 1978. The model results for CA 12 in the 
Spring suggested that in periods of rapid change (i.e., 
spring transition) a better estimate of the surface heat 
flux is required. The results for CA 1 emphasized this. 

It was found from the hindcast of CA 1 that the heat 
flux correction fields derived by Elsberry et al. (1982) 
were Peeisatisfactory for June to August 1976. A new 
surface heat flux correction field was derived specifically 
for 15 June to 18 August 1976. This new heat flux correction 
exactly accounts for the difference in the predicted (without 
correction) and analyzed temperature profile heat content 
relative to 200 m. With this more correct estimate of the 


surface heat flux, the model did produce an improved result 
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which was a fair representation of that analyzed for August 
1976. However, there were still errors in the model results. 
Perhaps some of the heat flux correction needs to be applied 
to the solar heat flux. There is also the possibility of 
model physics being the source of error. The Garwood model 
may benefit from the inclusion of entrainment shear production. 
This entrainment may reduce the sensitivity of the model by 
preventing overly shallow mixed layers from persisting 
during the spring and summer. 

The uncertainties in the heat flux correction fields 
are fairly large in the spring and summer. This limits 
tests of the ability of the model to predict anomalous 
Ocean conditions during the spring and summer. It is 
recommended that the model's capability for spring and 
summer predictions be further examined. 

The results of this study are generally encouraging. 
The most encouraging aspect is that a one-dimensional mixed 
layer model, such as the Garwood model, has the capability 
to predict anomalous conditions over long integration 
periods. The validity of a vertical mixing hypothesis for 
the formation and maintenance of cold anomalies in the 
North Pacific Ocean is also demonstrated. However, these 
predictions illustrate the crucial role of the surface heat 
and momentum fluxes in predicting the ocean mixed layer 


characteristics. The FNOC surface heat fluxes from the 
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hemispheric, primitive-equation model need to be improved, 
especially during the spring and summer. It is hoped that 
the new NOGAPS model at FNOC will provide these improved 


fields. 
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APPENDIX 


GARWOOD MODEL 

The Garwood (1977) model is a bulk or vertically 
integrated model for the oceanic planetary boundary layer 
(OPBL) or mixed layer. It is comprised of a closed system 
of seven equations. 


The entrainment buoyancy flux equation, 


—s 1/2 _ 
A <W > < Bie 
-bw(-h) = ———_;>————- (1) 


the horizontal component of turbulent kinetic energy, 


2 
3 eee = ERwin). AC | 
ego tv? ) = miu, 2AB 
2m ey 2) mee 
- m, (<B> - 3<w7>) (BE)? - [2 (ceo + 22 thy <B>; (2) 


il 


the vertical component of turbulent kinetic energy, 


i 3 @ — ae . 
oa (h <w >) = > hbw (-h) 5 ues, 
Wile 
ct m., (<E> = 3<w>) <E>i/2 
m m 
= (<E>l/2 + Beh) <E> ; 3) 
3 m4 
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the mean buoyancy and mean momentum equations 








a<B> a a 
h = = Bw(-h) - bw(o) + 2 9 (4) 
Pop 
RS “ = Gw(-h) - ew(0) + if<Coh (5) 


and the jump conditions at the bottom of the mixed layer 
(relating entrainment fluxes to the rate of deepending, 
and the changes in mean momentum and buoyancy at the base 


of the mixed layer), 


— _ dh 
cw (-~h) = yO and (6) 
— a dh 
bw (-h) = AB ae (a7) 


The time-dependent model forcing consists of the 
surface fluxes of momentum and buoyancy. In the present 
experiments, the contribution to the buoyancy due to 
salinity is neglected due to the lack of salinity data. 
Thus, the surface boundary conditions required to compute 
these fluxes are the total heat flux (Q.), the solar 
radiation (Q.) and the wind speed. Model outputs are the 
entrainment fluxes, turbulent kinetic energy, mixed layer 
depth and mixed layer temperature. 

The mechanism Garwood (1977) envisioned in the initial 


destabilization of the interface (between the mixed layer 
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and the denser water beneath) and the resulting entrainment 
Soma -OCal Kevin-Helmholtz instability (also known as a 
Pengemin (1963) class C instability). This instability is 
triggered by shear across the interface provided by the 

local turbulent eddies. The shear due to mean flow only 

is assumed to make a minor contribution to the achievement 

of this critical shear value; thus (5) can be neglected. 

By uSing the equation for the flux Richardson number evaluated 


at z = -h, 


Sema the jump conditions (6) and (7), one can find that 

shear production is a fixed fraction of buoyant damping in 
the entrainment zone. This zone may have a flux Richardson 
number greater than 1.0 and still possess enough turbulent 
kinetic energy for mixing to continue. In the Garwood model 
the shear production is only a secondary energy source for 
mixing, which can only be made available by entrainment 
initiated by another source. According to Garwood (1977), 
the most significant source of energy for mixing within an 
active entrainment zone is the convergence of flux of 


turbulent energy, 


3 E P 
— Gy tw Gig? 4 : 


ee, 





Thus, the critical parameter to determine the entrainment 
rate is not Rf but the ratio, P, of buoyancy flux to the 


convergence of energy flux: 


bw 


ho] 


3 
ne | WL +3) | 


The Garwood model uSes Separate vertical and horizontal 
equations (2) and (3) for the turbulent kinetic energy. 
This allows the convergence of turbulent kinetic energy 
flux to be included in the model. Mixed layer retreat 
occurs when the vertical component of turbulence is 
inadequate to transport heat, momentum and turbulence to 
the earlier depth of mixing. 

As explained by Tennekes and Lumley (1972), the dissi- 
pation rate can be estimated from the rate at which large 
scale eddies supply energy to the smaller scale eddies. 
This gives rise to a dissipative time scale 

ah 


ie = ase ; (10) 


M 


Where the vertically integrated averaged turbulent kinetic 


energy is defined as 


a ove <w > 
2 2 Z ; 
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ana the net rate of dissipation is 


<Ee> = f edz = if V =~ dz : (12) 


There are two time scales in the oceanic mixed layer. 
The first iS a convective time scale formed from the 
turbulent velocity and the length scale of the large scale 


Pempulent flow, t, = h/u, . This scale can be thought of 


i! 
as the time required for a large turbulent eddy to overturn. 


The second scale is imposed by planetary rotation, t. = l1/f 


ie 
In deeper boundary layers, planetary rotation turns the mean 


shear direction with depth and t, is the characteristic time 


2 
scale of the resulting vortex stretching. In Garwood's 
model, both of these time scales are incorporated into the 


parameterization of the dissipation time scale, which is 


defined as 


(13) 


The parameterization of the absorption of solar 
radiation used in the Garwood model is described by 
Gallacher et al. (1983) as a double exponential model for 
which one of the extinction depths approaches zero. 
Computationally, this reduces to absorbing a fixed fraction 


(l-r) of the flux of solar radiation, Qa in the first 
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meter, and exponential absorption for depths greater than 
ene) meter. 

Due to the process of Reynolds averaging in the turbu- 
lence equations and the subsequent formation of correlations 
(moments) of the fluctuating variables, there are more 
unknowns than equations. This problem can not be corrected 
at a fundamental level with additional equations since each 
new equation for a higher order of moment introduces addi- 
tional yet higher order moments. To close the system of 
equations, the chosen highest moments must be parameterized 
in terms of lower order moments and mean values. The Garwood 
model employs a second order closure scheme; thus the triple 
correlations (Second order moments) are parameterized in 
terms of zeroth order moments (mean values) and first order 
moments (autocorrelations and cross correlations). 

The process of parameterizing the high order moments 
generates constants of proportionality. In the set (1)-(3), 
my through mM, are constants of proportionality which must 
be determined from geophysical data. The values used in 
this study were m. = 210; P) = m,/™, = 1.0, P> = m.,/™, = 1.0, 
P3 = me /™4 =mOe5, 65 = 0.5 andy = 02001. The critical 
constants are M37 Par © and y. In order for the model to 
produce accurate forecasts it must be able to simulate the 
fall deepening, which is particularly sensitive to Mi the 
winter maximum mixed layer depth, which is sensitive to P3i 


the spring transition, which depends on the interaction of 
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M3 and the surface heat fluxes; and the summer maximum 
mixed layer temperature, which depends on r and y. The 
values of the critical constants listed above were 
empirically determined by Gallacher et al. (1983) to be as 
consistent as possible with the atmospheric forcing used 


imtnes study. 
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